Pressure-induced metal-insulator transition in LaMnC>3 is not of Mott-Hubbard type 
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Calculations employing the local density approximation combined with static and dynamical 
mean-field theories (LDA+17 and LDA+DMFT) indicate that the metal- insulator transition ob- 
served at 32 GPa in paramagnetic LaMnOa at room temperature is not a Mott-Hubbard transition, 
but is caused by orbital splitting of the majority-spin e g bands. For LaMnOs to be insulating at 
pressures below 32 GPa, both on-site Coulomb repulsion and Jahn- Teller distortion are needed. 

PACS numbers: 71.30.+h 71.20.-b 71.27. +a 
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Since the discovery of colossal magnetoresistance 
(CMR) 0, manganites have been intensively studied. 
The key to understand CMR is the high-temperature 
paramagnetic insulating-like phase, which is character- 
ized not only by an increase of resistivity with decreasing 
temperature, but also by unusual dynamical properties, 
such as low spectral weight at the Fermi level for a wide 
range of doping 0, S 0| • Theoretical understanding of 
this hole-doped paramagnetic phase remains incomplete, 
and CMR transition temperatures are lower than tech- 
nologically desirable. 

In this Letter, we shall focus on the pressure-induced 
insulator-metal (IM) transition in the undoped parent 
compound LaMnC>3 with configuration t\ g & g . This tran- 
sition occurs at room temperature, well above the mag- 
netic ordering temperature (Tjv=140K) , well below the 
cooperative Jahn- Teller (JT) temperature (T oo =740K at 
OGPa), and at a hydrostatic pressure of 32 GPa where 
the JT distortion appears to be completely suppressed 
The IM transition thus seems to be a bandwidth- 
driven Mott-Hubbard transition of the e g electrons and 
points to the dominating importance of the Coulomb re- 
pulsion between two e g electrons on the same site. This is 
supported by recent self-interaction-corrected local den- 
sity approximation (LDA) calculations, performed, how- 
ever, for the cubic structure and magnetically ordered 
states at low temperature . Structural distortions at 
K as functions of pressure were recently calculated with 
the LDA+C/ method 0- On the theoretical side, it has 
been an issue of long debate whether the JT distortion 
or the Coulomb repulsion is responsible for the insulat- 
ing behavior of LaMn03 at normal pressure. The high- 
pressure experiment H seems to favor the latter. 

Here, we shall study the room-temperature elec- 
tronic structure of LaMn03 at normal pressure and the 
jressure-induced IM transition by means of LDA + U 
g and LDA + dynamical mean field theory (DMFT) 
9] calculations. Upon going from the insulating to the 
metallic, high-pressure phase, we shall find that the or- 
bital polarization and the concomitant splitting of the 
two majority-spin e g bands are gradually reduced. The 
IM transition takes place when the bands start to overlap. 
Since this occurs within the (orbitally) symmetry-broken 



phase, this IM transition is not a Mott-Hubbard transi- 
tion. The Coulomb interaction, as well as the JT distor- 
tion are needed for a proper description of this transition 
and the insulating nature of LaMnOs. 

The orthorhombic crystal structure of LaMnOs at at- 
mospheric pressure is shown in Fig.^ The octahedra 
are elongated in the y direction (nearly parallel to b— a in 
Fig.GJ in subcells 1 and 3, and in the x direction (nearly 
parallel to b + a) in subcells 2 and 4. This JT distortion 
decreases linearly from 11% at OGPa, to 4% at 11 GPa, 
the highestpressure for which internal parameters were 
measured p| ■ The GdFe03-type distortion tilts the cor- 
ner sharing octahedra around the b axis and rotates them 
around the c axis, both in alternating directions. When 
the pressure increases from to 11 GPa, the tilt is re- 
duced from 12° to 8° and the rotation from 7° to 5°. 

The top right-hand side of Fig.[21 shows the paramag- 
netic LDA bands at normal pressure and the left-hand 
side those for a hypothetical, cubic structure with the 
same volume. The dashed blue bands were obtained 
with a lar ge b asis set of Ath-order muffin tin orbitals 
(NMTOs) Hi E3- Within the frame of the figure, we 
see the 3(x4) Mnt 2g and the 2(x4) Mne ff bands, while 
the O 2p bands are below. As a check, we also carried 




FIG. 1: (color). LaMnC>3 orthorhombic translational cell 
{Pbnm) and LDA-NMTO Mn e g crystal-field orbitals jl) (left) 
and |2) (right) of, respectively, lowest and highest energy. For 
the sake of clarity, the orbitals have been placed only in sub- 
cells 3 and 2; those in subcells 1 and 4 may be obtained by 
the LaO mirror plane perpendicular to the c (=z) axis. Since 
they have antibonding O 2p tails, orbitals |1) and |2) are di- 
rected, respectively, along and perpendicular to the longest 
Mn-0 bond. Red/blue indicates a positive/negative sign. 
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FIG. 2: (color). Top: Paramagnetic LDA band structure for 
orthorhombic (right) and hypothetical cubic (left) LaMnC>3 
at GPa plotted along the high-symmetry lines in the k x =k y 
plane. Energies are in eV, the k unit is 7r p"H|. and the k 
points marked are PMRXr in the cubic, and TYTZr in the 
orthorhombic BZ. The latter is folded in from the former and 
has the following smallest inequivalent reciprocal-lattice vec- 
tors: Q = 000, 110, 111, and 001. Dashed blue bands: large 
NMTO basis set; red bands: Mn e a NMTO basis set em- 
ployed in Eq. The N+l support energies Ei are shown 
at the right-hand sides. The zero of energy corresponds to 
configuration t\ g . Second row: and 11 GPa orthorhombic 
(OGPa cubic) e 9 bands folded out (in) to the (000,110)- 
BZ. The dimensionless band-shape parameter, X, is the LDA 
crystal-field splitting in units of the effective hopping integral 
t= \tdd<x\ ~W/6, both obtained from the NMTO Mn e 9 Wan- 
nier functions. Third row: As second row, but obtained by 
spin-polarized LDA+U for random spin orientations (room 
temperature). Bands are labelled by their main character. 
The zero of energy is the Fermi level. Bottom: Spectra cal- 
culated by LDA+DMFT. The full and dashed lines give the 
projections onto orbitals |1) and |2), respectively. 



out spin-polarized calculations for ferro- and antiferro- 
magnetic A-type orders and found full agreement with 
previous work [12j . Near the top of Fig.[5] — and contin- 
uing above it — is the La bd band which is pushed 2 eV 
up by pdo~ hybridization with oxygen when going from 
the cubic to the orthorhombic structure. This hybridiza- 
tion is a reason for the GdFeC>3-type distortion i but 
since it only involves Op orbitals perpendicular to the 



one which hybridizes with Mne 9 (Fig.^l, the La5d and 
Mne 9 bands hardly hybridize. Finally, the narrow band 
crossing the cubic Mn e g band is La 4/. 

The NMTO method can be used to generate minimal 
basis sets 0, El , such as the Mn e g basis whose orbitals 
are shown in Fig.^ This basis gives rise to the red solid 
bands in the topmost panel of Fig.|2 which are seen to 
follow the blue dashed bands exactly, except where the 
latter have avoided crossings with La bands. When their 
energy mesh is converged, the symmetrically orthonor- 
malized, minimal NMTO set is a set of Wannier func- 
tions. The e g NMTOs are localized by the requirement 
that a Mn e g orbital has no e g character on neighboring 
Mn atoms, and this confines the NMTO- Wannier func- 
tions to being essentially as localized as those in [T^j ]. 

Taking the Coulomb repulsion and Hund's exchange 
into account, three electrons localize in the ti g orbitals 
which we describe in the following by a classical spin S. 
These t2 g spins, which we assume to have random orien- 
tations at room temperature, couple to the e g electrons 
with strength 235=2.7 eV, as estimated by the splitting 
of the ej and bands in our ferromagnetic NMTO cal- 
culation (not shown). This results in the following low- 
energy Hamiltonian for the two e s -bands: 

ijmnaa' im 
+ U ^2 ^imt^iml + ^2( U ' ~ S t<t> -/)?W*W , (1) 



where hi 



is the LDA Hamiltonian in the represen- 



tation of the two (7n=l,2) e g NMTO- Wannier orbitals 
per site (Fig.^and two top rows of Fig|2J); u£ a , accounts 
for the rotation of the spin quantization axis (parallel to 
the t2 g spin) from Mn sites j to i. The second line de- 
scribes the Coulomb interactions between e g electrons in 
the same (17) and in different orbitals {U')\ J is the e g -e g 
Hund's rule exchange. We take £7=5 eV and J=0.75eV 
from the literature y|; by symmetry, U' = U — 2 J. These 
values are reasonable, also in comparison with those used 
for other transition-metal oxides. Whereas a larger U is 
appropriate when all five d degrees of freedom are treated 
[12| , the smaller value used for our e g Hamiltonian takes 
the screening by the t^ g electrons into account. 

We first solve {Q by DMFT using quantum Monte 
Carlo (QMC) simulations at room temperature. Previ- 
ous calculations for LaMn03 with electronic correlations, 
but simplified hopping integrals, include Ref. We ne- 
glected the (orbital) off-diagonal elements of the on-site 
Green function, forcing the DMFT density to have the 
same symmetry as the LDA crystal field (Fig.nj, a good 
approximation for e g systems. 

The spectral densities calculated using the observed 
crystal structures at and 11 GPa are shown at 
the bottom of Fig. [21 (green and blue, respectively). At 
GPa, we find strong orbital polarization, a 2 eV gap, 
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and, above 3eV, spectral densities which correspond to 
e gl and e g2 configurations antiparallel to the t2 3 spin at 
that site. Experiments for undoped LaMnC>3 show sim- 
ilar — albeit less sharp — gaps 0,0, 01 • At llGPa, we 
find the gap to be halved, but the orbital polarization to 
be just slightly reduced. 

Since for higher pressures the internal structural pa- 
rameters are unknown, we scaled the llGPa structure 
uniformly and calculated the LDA e g bandwidth W as a 
function of compression. The result is shown in blue on 
the left-hand side of Fig. OH On the right-hand side, we 
show the gap obtained from LDA+DMFT for various val- 
ues of W (blue solid line), keeping the LDA band shape 
constant. We see that the gap decreases, but does not 
close before W—7cV, which corresponds to a pressure 
way beyond the experimental 32 GPa. 

Let us now turn to the cubic phase (red), for which 
LDA+DMFT yields metallic behavior already at GPa 
(left-hand side of Fig.0) The spectral density ex- 

hibits a quasiparticle peak at £f, lower and upper Hub- 
bard bands at about +1.5 eV, and a peak due to the 
antiparallel-spin configurations around 4 eV. In the right- 
hand side of Fig.El we show how the quasiparticle weight 
Z decreases with volume expansion, and does not vanish 
before W=2.7 eV. One needs to pull the cubic crystal to 
make it insulating! Both types of structure, the cubic 
(X=0, red) and the one observed at 11 GPa (X=0.9i, 
blue) give grossly wrong pressures for the IM transition 
— way beyond the expected errors of the LDA and DMFT 
approximations. Hence, we can conclude that the IM 
transition occurs at a reduced (compared to the one at 
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FIG. 3: (color). Left: LDA e g bandwidths, W~6t, calculated 
as functions of compression, V/Vo- The top abscissa gives the 
experimental pressures, with 32 GPa marking the observed 
IM transition @. Right: 300 K LDA+DMFT (solid lines) 
and LDA+t/ (dashed lines) results calculated as functions of 
W. For insulators, we plot the DOS gap (top scale) and for 
metals the quasiparticle weights (bottom scale) . Full symbols 
indicate actual experimental structures and connecting black 
dotted lines are extrapolations. Each curve was calculated 
with fixed structure type (band shape): orthorhombic GPa 
(green, A=1.85), orthorhombic 11 GPa (blue, X=0.94), cubic 
plus crystal-field splitting (dark red, X=0.52), and cubic (red, 
X=0). 



11 GPa), but finite distortion [lfj. Such an intermediate 
distortion of X=0.52 (dark red line) indeed yields the IM 
transition at roughly the correct experimental pressure 
as seen by using the left hand-side of Fig-EJ to translate 
bandwidths into pressures. A finite distortion at the IM 
transition is also in accord with the observed noncubic 
crystal structure at the IM transition Besides the 
JT distortion, the Coulomb interaction is absolutely es- 
sential to describe insulating LaMn03. Without it, one 
would just obtain the metallic LDA bands shown in the 
second row of FigJ^l albeit with a width reduced by | 
due to the random spin orientation. 

A detaileded understanding of the band structure in 
the insulating phase may be obtained by treating the 
e g Hamiltonian Q in the static mean field approxima- 
tion (LDA+[7) as done by Ahn and Millis [l7|, but us- 
ing the accurate NMTO Hamiltonian. This is supported 
by the similarity of the resulting LDA+J7 gap-versus- 
W curves shown by the dashed green, blue, dark red, 
and red lines in the right-hand side of Fig.|21 to those 
obtained by DMFT. The LDA+C/ treatment overesti- 
mates the gaps 0] , but since it is easy to follow, we 
use it in the remainder of this Letter to explain — almost 
analytically — the development of the electronic struc- 
ture, starting with the cubic e g LDA bands, then turning 
on the orthorhombic distortions, and finally the spin and 
the e g -t2 g and e g -e g Coulomb interactions. 

In the cubic structure, the two e g bands never cross, as 
seen in the top left-hand side of Fig.El Since the domi- 
nating, indirect dpd, as well as the direct dd contributions 
to the hopping integrals (off-site elements of h im j n ) be- 
tween nearest neighbors in the z=0 plane are bonding 
(negative), the bottom of the band is at 000 and the top 
at 111. In the k x —k y plane the cigcnfunctions are the 
^3z 2 -i (k) and d x 2_ y 2 (k) Bloch functions, with the for- 
mer band lying lower. The hopping integrals are close 
to those of a first-nearest-neighbor tight-binding (TB) 
model with tdd8—0 and tddt?= — t- In this model, the cubic 
bandwidth is W=6t, and the band is symmetric around 
its center: £2 (k) = — £\ (k— [111] it). In the k x =k y =k 
plane, the dispersions of the d^-i and d x 2_ v 2 bands 
are, respectively, — t cos k — 2t cos k z and — 3icosfc. 

In order to understand the orthorhombic bands in the 
top right-hand side of Fig. [21 one should first note that 
mirroring an e g orbital in the LaO plane is equivalent 
to translating it along z. The 8 orthorhombic bands can 
therefore be folded out across the Brillouin zone (BZ) 
face fcs=f to the 4 green bands shown in the second 
row. The blue bands were calculated for the less dis- 
torted 11 GPa structure, and the 4 red bands are the 2 
cubic bands, folded into the Q= (110) it BZ 0. We see 
that the dominating effect of the distortion is to cou- 
ple the bands, and thereby to cut the e g band in two, 
which overlap by less than W. For k x =k y , there is no 
coupling between e x 2_ y 2 (k) and E3 Z 2_i (k) , but only be- 
tween £3 Z 2_ 1 (k) and e x 2_ y 2 (k — Q). This allows us to 
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fold the 4 bands out to merely 2, letting not only k z , 
but also k x =k y =k run from — n to n. In the TB model, 
these two bands are the eigenvalues of a 2 x 2 matrix with 
diagonal elements — t cos k — 2t cos k z + (A — \/3<5) /4 and 
3t cosfc — (A — V3£) /4 between, respectively, d 3z 2_ 1 (k) 
and d x 2_ y 2 (k — Q), and with the off-diagonal element 
- (v^A + S) /4. Here, tA/2 are the energies of the 
c?3j,2_ 1 and d x 2_ z 2 orbitals in cell 3 (Fig-QJ, —8/2 is the 

1 /2 

matrix element between them, and (A 2 + 5 2 ) =Xt is 
the crystal-field splitting. We see that the e g band is 
cut §t=jW above the bottom and |t below the top, 
with gaps of size Xt, so that the band overlap is merely 
(3-X)t + (X 2 ) . For the 11 and OGPa structures, 
respectively, (A, £)=(500, 110) and (849, 126) meV, and 
X=0.94 and 1.85. For large X, the density of states 
(DOS) gap is found to be {X-4)t + (X~ 2 ) . 

Inclusion of the spin and the e g -t2 g repulsion in the 
static mean field approximation splits the e g band by 
23S into e g and e g bands with spins locally parallel and 
antiparallel to that of the ti g spin. Because of the ran- 
dom orientation of the latter, the e g hopping inte gral s 
are independent of spin and reduced by the factor | |17fl . 
Since the coupling between the e g and e g bands is of or- 
der t' 2 /23S, the DOS gap remains (X' - 3.6) t' to order 
t'. Here, t'=M, X't'=Xt is the crystal-field splitting, and 
the value 3.6 for the IM transition was obtained numeri- 
cally; it lies between the small and large- A limits. 

Finally, the e gl -e g2 Coulomb repulsion, U" w U' — J = 
2.75 eV, splits the occupied and empty e g bands apart by 
the effective crystal field: 

(A 2 + PfJ* « (A 2 + <5 2 ) 1/2 + (n\ - n\) U" . (2) 

The orbital polarization, n\ — n\ = P (x' e ^J , is a 
band-structure function which increases linearly from 0, 
reaches 0.8 for X' e j:j=3.6, and saturates at 1 for large 
X' e j:f. As a consequence, Eq. written in the form: 
\X' e ff = X + P{X' e fj)U" /t is the self-consistency condi- 
tion for the effective crystal-field splitting and orbital po- 
larization as functions of the band-structure parameters 
t and X. This equation explains the dashed lines in Fig.|3 
and yields the band structures shown in the third line of 
Fig.0 Note that, in contrast to DMFT, the LBA+U ap- 
proximation gives a spontaneous orbital polarization of 
the cubic band structure at normal volume. 

In conclusion, the IM transition in LaMn03 at 32 GPa 
is not of Mott-Hubbard type. Rather, it is triggered by 
small distortions which create a crystal-field splitting, 
strongly enhanced by the Coulomb repulsion. For suf- 
ficiently large splitting, the majority-spin e g bands sepa- 
rate and LaMn03 becomes an insulator. Crucial are also 
the bandwidth reductions of | and %p arising from, re- 
spectively, the spatially uncorrelated directions of the t 2g 
spins at room temperature and the cutting by the effec- 
tive crystal field of the e g band in two subbands which 



overlap by less than the bandwidth. 
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Note added in proof. — An LDA+J7 study similar to 
ours, albeit for K and GPa, was reported in Ref. [IfSj . 



[1] R. von Helmolt et al, Phys. Rev. Lett. 71, 2331 (1993); 

S. Jin et al, Science 264, 413 (1994); P. Schiffer et al, 

Phys. Rev. Lett. 75, 3336 (1995). 
[2] A. E. Bocquet et al, Phys. Rev. B 46, 3771 (1992); M. 

Abbate et al, Phys. Rev. B 46, 4511 (1992); A. Chainani 

et al, Phys. Rev. B 47, 15397 (1993); T. Saitoh et al, 

Phys. Rev. B 56, 8836 (1997). 
[3] J.-H. Park et al., Phys. Rev. Lett. 76, 4215 (1996). 
[4] Y. Okimoto et al, Phys. Rev. Lett. 75, 109 (1995); K. 

Takenaka et al, J. Phys. Soc. Jpn. 68, 1828 (1999); K. 

Tobe et al, Phys. Rev. B 64, 184421 (R) (2001); M. A. 

Quijada et al, Phys. Rev. B 64, 224426 (2001); N. N. 

Kovaleva et al, Phys. Rev. Lett. 93, 147204 (2004). 
[5] I. Loa et al, Phys. Rev. Lett. 87, 125501 (2001). 
[6] G. Banach and W. M. Temmerman, J. Phys.: Condens. 

Matter 16, S5633 (2004); H. Zenia et al, New J. Phys. 

7, 257 (2005). 

[7] G. Trimarchi and N. Binggeli, Phys. Rev. B 71, 035101 
(2005). 

[8] V. I. Anisimov et al, Phys. Rev. B 44, 943 (1991). 

[9] V.I. Anisimov et al, J. Phys.: Condens. Matter 9, 7359 

(1997); A.I. Lichtenstein and M.I. Katsnelson, Phys. Rev. 

B 57, 6884 (1998); K. Held et al, Psi-k Newsletter 56, 5 

(2003). 

[10] O. K. Andersen and T. Saha-Dasgupta, Phys. Rev. B 62, 
R16219 (2000); E. Zurek et al, ChemPhysChem 6, 1934 
(2005). 

[11] E. Pavarini et al, New J. Phys. 7, 188 (2005); E. Pavarini 
et al, Phys. Rev. Lett. 92, 176403 (2004). 

[12] S. Satpathy et al., Phys. Rev. Lett. 76, 960 (1996); W. E. 
Pickett and D. J. Singh, Phys. Rev. B 53, 1146 (1996); I. 
Solovyev et al, Phys. Rev. Lett. 76, 4825 (1996); Phys. 
Rev. B 53, 7158 (1996); H. Sawada et al, Phys. B 56, 
12154 (1997); I. S. Elfimov et al., Phys. Rev. Lett. 82, 
4264 (1999). 

[13] N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 
(1997). 

[14] W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989); A. Georges et al, Rev. Mod. Phys. 68, 13 (1996). 

[15] Y. Motome and M. Imada, J. Phys. Soc. Jpn. 68, 16 
(1999); K. Held and D. Vollhardt, Phys. Rev. Lett. 84, 
5168 (2000); V. Ferrari and M. J. Rozenberg, Mod. Phys. 
Lett. B 15, 1031 (2001); Th. Pruschke and M.B. Zolfl, 
Adv. in Solid State Phys. 40, 251 (2001); M. S. Laad et 
al., New J. Phys. 6, 157 (2004). 

[16] LDA+DMFT tests for cubic LaMn0 3 at GPa and 600 K 
gave no spontaneous orbital polarization. 

[17] K.H. Ahn and A.J. Millis, Phys. Rev. B 61, 13545 (2000). 

[18] G. Sangiovanni et al, cond-mat/05 11442 

[19] W.-G. Yin et al, Phys. Rev. Lett. 96, 116405 (2006). 



